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Abstract 

We have studied quasinormal modes of spinor massless Dirac field perturbations 
of a black hole in massive gravity. The parameters of the theory, such as the mass of 
the black hole, the scalar charge of the black hole, mode number and the multipole 
number are varied to observe how the corresponding quasinormal frequencies change. 
We have also used the Pdschl-Teller approximation to reach analytical values for 
the frequencies of quasinormal modes for comparison with the numerically obtained 
values. Comparisons are done with the frequencies of the Schwarzschild black hole. 
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1 Introduction 

Massive gravity is an alternative theory to General Relativity where the graviton, 
which is a spin two held, acquire a mass. There are several reasons to study massive 
gravity theories: one of the reasons being able to explain the acceleration of the uni¬ 
verse without introducing the cosmological constant or dark energy. It is speculated 
that by introducing mass for the graviton, the gravity can be modihed at the infrared 
in a such a way as to produce the acceleration of the universe. 

The hrst attempt to include mass for the graviton was done by Fierz and Pauli in 
1939 [T]. In comparison with GR where the graviton has two degrees of freedom, in 
a massive gravity theory the graviton has hve degrees of freedom. 

During the last few years, signihcant progress has been done on massive gravity 
theories. For example, a particular massive gravity theory called dRGT theory 00 is 
found to be free from Boulware-Deser ghosts. Two other well known massive gravity 
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theories which are also know to be free from ghosts are, DGP model [1] and the “new 
massive gravity theory” in three dimensions [5]. There are many interesting works 
related to massive gravity in the literature and we find it difficult to discuss all; we will 
direct the reader to excellent reviews on the topic by de Rham [6] and Hinterbichler 
[7] instead. 

In this paper, we focus on a massive gravity theory which is a Lorentz violating 
theory. In this model, Lorentz symmetry is broken spontaneously by four scalar 
fields called Goldstone bosons. These field are coupled to gravity through derivative 
coupling. When Lorentz symmetry is broken spontaneously, the graviton acquire a 
mass very similar to the Higgs mechanism. A nice review of Lorentz violating massive 
gravity theory can be found in [8] [9]. 

The Lorentz violating theory of massive gravity considered in this paper is de¬ 
scribed by the following action: 


S= 


R 

.Ibvr 




( 1 ) 


Here R is the scalar curvature of the space-time geometry and is a function of 
X and hL*L The functions X and are functions of scalar fields and are 

defined as. 


“ A4 

A4 A^x 


( 2 ) 

(3) 


The theory described by action in eq.([T]) is the low-energy effective theory below the 
ultra-violet cutoff A and a perturbative analysis on the theory calculate the value of A 
to be in the order of ^mMpi where m is the graviton mass and Mpi the Plank mass [10] 
[8][n][i2]. The scalar helds are responsible for spontaneously breaking Lorentz 

symmetry. When symmetry is broken, the scalar fields acquire a vacuum 

expectation value and are called Goldstone fields. 

Studies of perturbations and quasi normal modes (QNM) of black holes have a 
rich history. The founders of black hole perturbations were Regee and Wheeler where 
they established the equations of axial perturbations of the Schwarzschild black hole 
in 1957 [I3]. Thereafter, the contributions from Zerilli, Vishveshwara, Press and 
Teukolsky put a hrm foundation on the area of research on black hole perturbations 
as an important held of research in gravitational physics. A nice review on black hole 
perturbations and current methods used can be found in [H] . QNM are a set of modes 
arising during the intermediate stage of a black hole perturbation. The frequencies 
of QNM are complex and only depend on the properties of the black hole such as the 
mass, spin and the charge. Studies of QNM frequencies has attracted great attention 
due to variety reasons: one of the prominent being the current attempts to observe 
gravitational waves. If current gravitational wave detectors can detect signals emitted 
from black holes, then we will be able to specify their properties of such objects in the 
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universe [15]. An exhaustive reasons as to why we need to study QNM frequencies 
are given in the nice review by Konoplya and Zhidenko [T6] . 

In this paper we focus on the spin massless Dirac spinor held perturbation 
of a black hole in massive gravity. There are several works which have focused on 
the QNM of Dirac perturbations of black holes: Dirac perturbations of the Lifshitz 
black hole in 4 dimensions were studied in im. Scalar and Dirac perturbations of 
the Kerr-Newmann-de Siter black hole were analyzed by Konoplya in [18]. Spinor 
perturbations of regular black holes where nonlinear electrodynamics are coupled to 
gravity were studied by Li and Ma in [19]. Decay of spinor half held around a Born- 
Infeld black hole was studied by Fernando poi. 

The paper is organized as follows: in section 2, the black hole in massive gravity 
is introduced. In section 3, the equations for the massless Dirac held of spin ^ is 
developed and presented. In section 4, the WKB approach is employed to compute 
QNM frequencies by varying the parameters in the theory. The Pdschl-Teller approx¬ 
imation is used to hnd analytical values for uj in section 5. Finally in section 6 the 
conclusion is given. 


2 Black holes in massive gravity 

In this section we will present the basic characteristics of the black hole in massive 
gravity considered in this paper. A detailed derivation can be found in [21] and [22] . 
The geometry is described by the metric given by, 

ds^ = —fir)dt^ -I- — -I- r'^idO"^ -I- sin^Odcf)^) (4) 

f{r) 


where, 


/(O 


The scalar helds are given by. 


2M 

r 


Q 


$0 = A^t + z{r )); <l>* = AV 


(5) 

( 6 ) 


where 


z(r 


= ± 


dr 

fir) 


1 - fir) 


QA(A-l) 1 

12m2 r^+2 



(7) 


Here m is the mass of the graviton and A is a positive constant. A has to be greater 
than 1 in order for the solutions to be asymptotically hat and M to be the ADM 
mass. Q is a scalar charge and represents a gravity theory with a massive graviton. 
The function jr for this particular black hole solution is given by. 


jr 


12 

T 


i+p. 


(pI - 3piP2 - Qpi + 2p3 - 12 ) 


( 8 ) 
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where, 


Pn = Tr{W^) (9) 

There are two possibilities for Q when M > 0. When Q > 0, the geometry is very 
similar to the Schwarzschild black hole with a single horizon. The function /(r) for 
this case is given in Fig.([T]). The horizon radius for the black hole in massive gravity 
is larger than the one for the Schwarzschild black hole. When A —>■ cxo, r/j —)• 2M 
which is the value for the Schwarzschild black hole. 

When Q < 0, the geometry is similar to the well known Reissner-Nordstrom 
charged black hole. There could be two, one or no horizons depending on the param¬ 
eter of the theory. The function /(r) for this case is given in Fig.(|2]). 

When the mass of the black hole M is. 


M, 


critical — 


A|Q| 


l/A 


A- 1 


A-1 

A 


( 10 ) 


the horizons merge. For M > McriUcai, there will be two horizons. For M < Mcriucai, 
there wont be any horizons and there will be a naked singularity. 



2 4 6 8 10 


r 

Figure 1: The hgure shows /(r) vs r. Here M = 1 and Q = 1. 
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Figure 2: The figure shows /(r) vs r. Here Q = —2 and A = 2. 


The Hawking temperature of the black hole is given by, 



df{r) 


dr 


r=rh 


^ (‘m Q\ \ 
471 \rl ^ r^+V 


( 11 ) 


The temperature is plotted varying the mass M and Q in Fig.(j3]) and Fig.(j4]). When 
M increases, the temperature decreases. This behavior is similar to the behavior of the 
temperature of the Schwarzschild black hole. When Q is increased, the temperature 
increases to a maximum and then decreases. Thermodynamics and phase structure 
of the massive gravity black hole analyzed in this paper was addressed by Capela and 
Nardini [23] . 



Figure 3: The figure shows Temperature vs M. Here Q = 1 and A = 3. 
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Figure 4: The figure shows Temperature vs Q. Here M = 1 and A = 3. 


3 A massless, spin | Dirac field around a black 
hole 


In this section, the Dirac equation around the black hole considered in this paper is 
presented. The equation for a Dirac spinor field of spin A in a curved back ground is 
given by, 

re:{d, - r,)e = o (12) 

Here, e“ is the tetrads and is the inverse of the tetrads. The metric tensor ga/s of 
the space-time considered and the tetrads e“ are related by, 

9a(3 (^ 3 ) 

Here, rjab = (—1,1,1,1) is the metric of the fiat space. The 7 “ matrices appearing in 
eq. ffT 2 |) are defined by. 


(^ 0 y 0 

0 -i y ’ ' 0 J 


Here, cr“ are the well known Pauli matrices given by. 


a 


1 


1 

0 



a 


0 -i\ 

i 0 y ’ 


a 



The 7 “ matrices satisfy the anti-commuting relations. 




( 14 ) 


( 15 ) 


( 16 ) 


The spin connections T^ in eq. flT^ are defined as. 
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r„ = -jl7“.7'Kei„ 


with, 




eb/3 


are the Christoffel symbols. 


For the metric given in this paper, the tetrads are given by. 


e“ = diag r, rsinOj 


( 17 ) 

(18) 

(19) 


Therefore, the spin connections for this space-time is given by, 

P /' 0 1 

Tt = —7 7 
= 0 

P v7 1 2 
Te = ^7 7 

= ^7^7^ + ^sm07S3 

To facilitate computations, the function ^ can be redehned as, 

$ 




yi/i 


Then the Dirac equation flT^ simplihes to. 


7 O 

V7 dt 




—+ - 
dr r 


$ + 


7^ 


' d cot9 ' 


<h + 


70 




rsinO dip 


= 0 


( 20 ) 


( 21 ) 


( 22 ) 


Before proceeding further, some clarihcations are required. Since the field considered 
here is massless and spin |, the allowed solutions to the Dirac equation are circularly 
polarized. More details of this argument can be found in the paper by Brill and 
Wheeler [2l]. According to Brill and Wheeler [2l], the spinors considered have right 
handed circular polarization and the allowable spin states satisfy the identity. 


Here 75 is given by. 


(1 - ^ 75 )$ = 0 


= 7*^7^7^7^ 


(23) 

(24) 


The eq. fl23l) gives a simplihed set of components for $ as. 


$ = 


^2{t,r,e,ip) 
^i{t,r,e,ip) 

\ ^2{t,r,e,ip) ) 


( 25 ) 
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Hence the Dirac equation yield two identical set of equations each coupling <hi and 
^ 2 - The components $1 and $2 are redehned as follows: 


) \^a2ie,ip) ) 

Substituting eq. fl^ to the Dirac equation ([22]) yields, 


zur 

yf 


K 






14 + 


' cot9 


dr J G 


-02 + 


sin9 dcp 


da2 da2 

+ 



(26) 


(27) 


'iur r:dG\ 1 / cot9 


doi doi 

+ 


sin9 dip 89 ) 02 


= 0 


(28) 


The angular part of the above equation can be solved in terms of spin-weighted 
spherical harmonics, sY)m- A lengthy description of spin weighted spherical harmonics 
can be found in |25] and [26]. One can dehne two operators 9+ and 8- as. 


d+ 


cot9 i 8 5 \ 

2 sin9 8ip 89 ) 


(29) 


cot9 i 8 8 \ 

2 sin9 8<p 89 ) 

It was shown in [25] [26] that the above operators act on spin weighted spherical har¬ 
monics .jT/m as ladder operators. In particular if the spin s = ^, the above operators 
gives the following relations; 




8+{_lYiJ = {l + l-)iYim 

2 ^ ^ 

(31) 

d-UYlm) = -{l+h-lYlm 

2 2i 

(32) 


Note that / < |s| and I has to be half integer. Hence I can be written in terms of p 
which is a positive integer as, / = p — Here p will be called the multipole number 
in the rest of the paper which can take any positive integer. 

If the functions ai{9,ip) and 0 : 2 ( 6 *, (p) are chosen as spin weighted spherical har¬ 
monics given by, 

«i = (33) 

02 = -iYim (34) 

the eq. (l27)) and eq. fl28]) can be simplihed to be. 


dr* 


iuG + W{r)K = 0 


(35) 













( 36 ) 


— + iojK + W{r)G = 0 
dr* 

r* is the “tortoise” coordinate given by, 

dr* = y (37) 

and the fnnction W{r) is, 

W{r) = (38) 

r 

Two new fnnctions (5± can be dehned as, 

P± = K±G (39) 

and eq. fl35|l and eq. fl36|) can be deconpled as, 

^ + (n;2 - = 0 (40) 

Here, are related to W (r) as, 

= =F/ (41) 

In a paper by Anderson and Price |2^ it was discnssed that two potentials and 

Vbirac given above which are related, will prodnce the same physical consequences. 
Hence, both potentials will produce the same QNM spectra. Therefore, we will use 
only VDjj-ac of our computations in the rest of the paper and will be referred to 

as just VDirac(r). 

The effective potential VDiradj") depends on four parameters: M, Q, fi and A. In 
Fig. ([5]), VDirac{r) is plotted as a function of r by varying A. When A increases, the 
height of the potential increases. In Fig. ([ 6 ]), Voiradr) is plotted as a function of r 
by varying /r. When /x increases, the height of the potential increases. In Fig. ([7]) 
Voiracd) is plotted as a function of r by varying M. When mass increases, the height 
of the potential decreases. In Fig. ([8]) , Voiracd) is plotted as a function of r by 
varying Q. When Q increases, the potential height decreases. However, the potential 
for the Schwarzschild black hole (with Q = 0) is higher than any one of them. The 
behavior of the potential Vmrac is similar to the behavior of the effective potential of 
the massless scalar held around the massive gravity black hole [37] . 
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Figure 5: The figure shows VDiracir) vs r. Here M = 2,/i = 3 and Q = 1. 



Figure 6: The figure shows Voiracij') vs r. Here M = 1, A = 1 and Q = 1. 
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Figure 7: The figure shows Voiracir) vs r. Here A = 1,/i = 2 and Q = 1. 



Figure 8: The figure shows Vrurad.'^') vs r. Here M = 0.8, A = 1 and fi = 2. 


4 Computation of QNM frequencies 

QNM of a black hole for a particular perturbation are given by the solutions to the 
corresponding wave equations. In this paper, the QNM are given by the solution to 
the Dirac wave equation in eq lHnjl . To seek solutions boundary conditions has to be 
imposed; for asymptotically fiat black hole, such as the one considered in this paper, 
the boundary conditions are, ingoing waves at the horizon and out going waves at the 
spatial infinity. In this work we are particularly interested in frequencies of QNM. 
QNM frequencies are complex and do not form a complete set. For stable black holes, 
the imaginary part of the frequency u is expected to be negative. The modes of the 
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perturbations are labeled by an integer n. The fundamental mode is described by 
n = 0 mode. Since the largest value of uj is for n = 0, we will mostly be interested 
on ooi with n = 0. 

Analytical solutions to the wave equations resulting from black hole perturbations 
are rare. There are few examples we would like to mention here. QNM frequencies of 
dilaton black holes in 2+1 dimensions have been found by Fernando in [28] [29] [30] . 
Some other examples with exact QNM frequencies are in [31] [32] [33]. 

Due to the inability to solve the wave equations and frequencies analytically, there 
are numerous numerical methods developed to compute QNM frequencies. The reader 
is referred to the review by Konoplya and Zhidenko [I6] as an excellent source of ref¬ 
erence for various methods. Out of the methods used in literature, the WKB analysis 
is one of the most popular methods to hnd QNM frequencies. WKB approach was de¬ 
veloped by Iyer and Will [3l] and subsequently developed to sixth order by Konoplya 
[35] . This method has been employed to compute QNM frequencies corresponding to 
black hole perturbations in several cases. Examples are found in [20] [30] [3Z] • In this 
method, the QNM frequencies are given by the expression, 

= ~‘i\J ^ 2 ) ^i'^rnax) (42) 

Here, r^ax is the place where V (r) is maximum and V''{r) is the second derivative of 
the potential. Expressions for Dj can be found in [35] . 

In this paper, we used the WKB approximation described above to compute the 
QNM frequencies. We used all 6 orders for the computation. WKB approximation is 
accurate when the multipole number fi is greater than the overtone number n [35] . 

We have computed u for the Dirac held by varying the parameters in the theory, 
M, A, Q, n and fi. All uji values were negative; therefore, the black hole is stable under 
Dirac held perturbations. We have presented results in tables as well as hgures. Notice 
that uj is negative and we have only written the value without the negative sign in 
the table. We have also used only the value in plotting the graphs. 

First, let us discuss the ehects of A on a;. The values are given in Table 1. For 
small values of A, both ur and uj increases with A. When A is increased further, ojr 
approaches a constant value. The behavior is similar for both u = 0 and u = 1, and is 
presented in Fig. ([9]). For cu/, increasing A further increase it to a maximum and then 
decreases to a constant value. This behavior is demonstrated in Fig. fITU]) . The stable 
values reached by both ojr and uj corresponds to the values for the Schwarzschild 
black hole given by, uisch = 0.09133— i0.04747(/i = 1) and 0.1779 —'i0.1487(/i = 2). In 
general, the Dirac held decay slower around the massive gravity black hole for small 
A. There is a range of A for which the Dirac held decay faster than the Schwarzschild 
black hole. Hence the signihcance of A on the decay rates are only for a short range as 
is obvious from the graphs. This behavior is not surprising since the metric function 
/(r) 1 — ^ for large A values. The behavior of the scalar held around this black 

hole was similar as shown in [37] . 
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When ojR and oji were studied by varying the mass of the black hole as given in 
Table 2, both decreases with increasing mass M as shown in Fig. lfTTD and Fig. ([12]). 
Hence, small black holes are more stable compared to large black holes. The scalar 
held behaves similarly 133. 

When Q parameter is varied, ur decreases with Q. On the other hand uj increases 
to a mamximum and then decreases as shown in Table 3 and Fig. flT^ . Hence there 
is a maximum value of Q which gives more stability. Interestingly, the scalar held 
demonstrated diherent behavior by increasing uj when Q is increased |37] . 

In Fig. llTTll and Fig. llTHll . ur and uj are plotted against the multipole number /i. 
The values are given in Table 4. When /i is increased, ur increases linearly. For 
both n = 0 and n = 1, this behavior persists. On the other hand, ojj decreases to a 
constant value when I is increased. The behavior is similar for n = 0 and n = 1. 

When ojr and ujj is plotted by changing the mode number n, ur decreases with n 
and ujj increases with n. Also, uj has a linear relation with n. This is demonstrated 
in Fig. llTB]) and Table 5. 

It was observed that the temperature of the black hole depend on uj linearly 
for asymptotically anti de Sitter black holes |38]. Hence, we also plotted ui vs the 
temperature of the black hole by varying Q. A linear relation seems to exist between 
the temperature and uji 


A 

Ur (n =0, /i = 1) 

Ui{ n =0, /i = 1) 

cu/j ( n =1, /i = 2) 

Uj ( n =1, fi = 2) 

2 

0.08759 

0.04663 

0.1700 

0.1468 

2.2 

0.08864 

0.04708 

0.1721 

0.1481 

2.4 

0.08940 

0.04738 

0.1736 

0.1488 

2.6 

0.08994 

0.04757 

0.1747 

0.1492 

2.8 

0.09032 

0.04768 

0.1756 

0.1494 

3 

0.09060 

0.04774 

0.1762 

0.1494 

3.2 

0.09079 

0.04777 

0.1767 

0.1494 

3.4 

0.09092 

0.04777 

0.1770 

0.1493 

3.6 

0.09102 

0.04775 

0.1773 

0.1492 

3.8 

0.09109 

0.04772 

0.1775 

0.1491 

4 

0.09115 

0.04769 

0.1776 

0.1490 

4.2 

0.09119 

0.04766 

0.1777 

0.1490 

4.4 

0.09121 

0.04762 

0.1778 

0.1489 

4.6 

0.09124 

0.04759 

0.1778 

0.1488 

4.8 

0.09127 

0.04756 

0.1779 

0.1488 

5.0 

0.09128 

0.04754 

0.1779 

0.1487 

5.2 

0.09130 

0.04752 

0.1779 

0.1487 

5.4 

0.09131 

0.04750 

0.1779 

0.1487 

5.6 

0.09132 

0.04748 

0.1779 

0.1487 

5.8 

0.09132 

0.04747 

0.1779 

0.1487 

6.0 

0.09133 

0.04747 

0.1779 

0.1487 
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Table 1: QNM frequencies for various A values. Here M = 2, Q = 1. 


M 

Ur (n =0, /i = 1) 

Uj{ n = 0 , /i = 1 ) 

Ur {n =1, IJ, = 2) 

Uj ( n = 1 , fi = 2) 

1 

0.1735 

0.09679 

0.3321 

0.3061 

1.5 

0.1196 

0.06406 

0.2320 

0.2005 

2 

0.09060 

0.04774 

0.1762 

0.1494 

2.5 

0.07276 

0.03809 

0.1416 

0.1193 

3 

0.06074 

0.03170 

0.1183 

0.0993 

3.5 

0.05211 

0.02716 

0.1015 

0.0850 

4 

0.04562 

0.02375 

0.0889 

0.0744 

4.5 

0.04056 

0.02111 

0.0790 

0.0661 

5 

0.03651 

0.01900 

0.0711 

0.0595 

5.5 

0.03320 

0.01727 

0.0647 

0.0541 

6 

0.03043 

0.01583 

0.0593 

0.0496 


Table 2: QNM frequencies for various M values. Here A = 3, Q = 1. 


Q 


Ui 

Temperature 

0 

0.09132 

0.04747 

0.0198944 

1 

0.09060 

0.04774 

0.0201830 

2 

0.08991 

0.04797 

0.0204320 

3 

0.08928 

0.04815 

0.0206480 

4 

0.08869 

0.04828 

0.0208380 

5 

0.08815 

0.04837 

0.0210050 

6 

0.08764 

0.04841 

0.0211536 

7 

0.08717 

0.04842 

0.0212854 


Table 3: QNM frequencies for Q values. Here M = 2, A = 3 and /i = 1. 


fi 

Ur (n =0) 

Ui{ n = 0 ) 

Ur ( n =1) 

ui ( n = 1 ) 

2 

0.1826 

0.04750 

0.1700 

0.1468 

3 

0.2760 

0.04746 

0.2671 

0.1444 

4 

0.3689 

0.04744 

0.3621 

0.1434 

5 

0.4617 

0.04743 

0.4562 

0.1430 

6 

0.5544 

0.04743 

0.5498 

0.1426 

7 

0.6471 

0.04742 

0.6431 

0.1426 

8 

0.7397 

0.04742 

0.7362 

0.1426 

9 

0.8323 

0.04742 

0.8292 

0.1425 

10 

0.9249 

0.04742 

0.9221 

0.1424 

11 

1.0175 

0.04742 

1.0150 

0.1423 

21 

1.9431 

0.04742 

1.9418 

0.1423 

31 

2.8686 

0.04742 

2.8677 

0.1423 
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Table 4: QNM frequencies for various /i values. Here M = 2, A = 2 and Q = 1. 


n 


UJi 

0 

1.01748 

0.04742 

1 

1.01496 

0.14241 

2 

1.00996 

0.23784 

3 

1.00256 

0.033402 

4 

0.099287 

0.43122 

5 

0.98104 

0.52972 

6 

0.96727 

0.62975 


Table 5: QNM frequencies for various overtone mode number n values. Here 
M = 2, A = 2, /i = 11 and Q = 1. 




Figure 9: The figure shows ujr vs A. Here M = 2,Q = 1 and /x = 1, 2. 




Figure 10: The hgure shows uj vs A. Here M = 2,Q = 1 and /x = 1, 2. 
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Figure 11: The figure shows ojr vs M. Here A = 3, Q = 1 and /i = 1, 2. 




M M 


Figure 12: The figure shows oji vs M. Here A = 3, Q = 1 and /x = 1, 2. 
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Figure 13: The figure shows ojr and uj vs Q. Here M = 2, A = 3 and /i = 1. 
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Figure 14: The figure shows ojr vs /i. Here M = 2, Q = 1 and A = 2. 



Figure 15: The figure shows oji vs fi. Here M = 2, Q = 1 and A = 2. 
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Figure 16: The figure shows ojr and oji vs n. Here M = 2, Q = 1, /x = 11 and A = 2. 
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Ul 


Figure 17; The figure shows Temp vs ui. Here M = 2,/i = 1 and A = 3. Here Q is 
varied 


5 Analytical values of uj with Poschl-Teller ap¬ 
proximation 

One can obtain analytical values of u by approximating the effective potential with the 
well known Poschl-Teller potential. Ferrari and Mashhoon [39], applied this method to 
obtain exact results for the Schwarzschild black hole, Reissner-Nordstrom black hole 
and the Kerr black hole. Since the wave equation with the Poschl-Teller potential 
can be solved exactly, one can find analytical formulas for the QNM frequencies u. 
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Here, we will apply this method to obtain approximate values for QNM frequencies 
of the Dirac held with large A. 

In the the Pdschl-Teller approximation, the effective potential is approximated by. 


V 


Vo 

Cosh'^/3(r^ — r^o) 


(43) 


Here r* is the tortoise coordinate described in eq. fl57|) . 
potential has its the maximum: hence ^ = 0 at r* = r*o. 
given by. 


/3^- 

Vo 


1 (fV 

2Vo drl 
= 


r=r*o 


r^o 


) 


r*o is the point where the 
The variable /3 and Vq are 

(44) 

(45) 


It was shown in [?] [39] that the QNM frequencies u for the above potential is given 

by, 

a; = ±^Vo - )574 - i/3(n + 1/2) (46) 

Due to the complicated nature of the potential in eq. flTT]) . we will study the QNM 
frequencies at the eikonal limit (large /i). Then, the dominant term in the effective 
potential is the one proportional to jj?. Hence, 


V{r) ^ 



(47) 


Since we are studying the wave equation for large A, the function /(r) can be approx¬ 
imated as, 

fir) Ri 1 — — Qk (48) 

r 

Here = j and Q is chosen such that / will be dimensionally correct. When A —)■ 
cx), / ~ 1 — —. Hence for large A, the horizon gets closer to the horizon of the 
Schwarzschild black hole which is what is expected as given in Fig.([T|). 

Now, = (^) (Jr) = 0 lead to the solutions. 


3M 

(1 - Qk) 


The hnal results for /? and W are given as. 


(49) 


2 (1 - 


/S = 

Vo = ^^^ 


27M2 

(1 - Qk)^ 

27M2 


(50) 

(51) 
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Hence u is given as, 


{1-Qk)‘^ IJ? 1 ^{I-Qk^ 

~lh/3M~ \ (1 -Qk) ~4~ ' 3\/3M 


(52) 


It is clear that when A —?• cxd(k 0), a; reaches the valne for the Schwarzschild black 
hole obtained by Ferrari and Mashhoon in |39j . 

By observing the u given above, one can explain the behavior of u obtained with 
the WKB approach. For example, when /x becomes large, the above approximation 
gives ujji oc /X which was observed in Fig. ffn|) . On the other hand ox/ is independent 
of /i which is what was demonstrated in Fig. ffTSjl . When /t —)■ 0(A —)■ oo), uj —)■ 
(^Schwarzschild- This is dear from Fig. ([9]) and Fig. flTOl) . When the mass M is increased, 
both ojr and ox/ decreases from eq. flS^ . This is the same behavior observed in Fig. ffTTl) 
and Fig. ffT^ . 


6 Conclusion 

We have studied QNM frequencies of a spin massless Dirac held perturbations of 
a black hole in massive gravity. First we computed frequencies ax with the 6th order 
WKB approximation. The parameters in the theory. A, M, Q and /x are varied to 
obtain the relations with ax. For small values of A, both ur and ox/ increases. When 
A gets larger, ojr reaches a stable value; ox/ keep increasing to a maximum and then 
decreases to a stable value. The stable value corresponds to the ox of the Schwarzschild 
black hole. Similar behavior was observed for the massless scalar held perturbations 
around the black hole |37]. 

When the mass M is increased, both ur and ox/ decreases. Hence, small black 
holes are more stable. When Q increases, ojr decreased. On the other hand ox/ 
increased to a maximum before falling oh for large Q. Hence there is a maximum 
value of Q for which the black hole is most stable. 

When the multipole number /x is increased, oxr increase linearly; ox/ decreases and 
reach a stable value. This behavior is consistent for the mode values rx = 0 and n = 1. 
When the mode number n is increased, ujr decreases; ox/ increases linearly with n. 

We used Pdschl-Teller approximation to hnd an analytical values for ox when A 
is large. Our computation was done for large /x. The analytical ox obtained in this 
method verihed the numerical values obtained with the WKB approach. 

As future work, it would be interesting to hnd the stability under electromagnetic 
perturbations. Also, the cross sections of helds of spin 0, 1 would be an interesting 

aspect to study. 

Acknowledgments: The author wish to thank R. A. Konoplya for providing the 
Mathematica hie for WKB approximation. 
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